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Abstract 

A new simplified approach for teaching electromagnetism is presented using the formahsm of 
geometric algebra (GA) which does not require vector calculus or tensor index notation, thus pro- 
ducing a much more accessible presentation for students. The four-dimensional spacctime proposed 
is completely symmetrical between the space and time dimensions, thus fulfilling Minkowski's orig- 
inal vision. In order to improve student reception we also focus on forces and the conservation of 
energy and momentum, which take a very simple form in GA, so that students can easily build 
on established intuitions in using these laws developed from studying Newtonian mechanics. The 
potential formulation is also integrated into the presentation that allows an alternate sohition path, 
as well as an introduction to the Lagrangian approach. Several problems are solved throughout 
the text to make the implementation clear. We extend previous treatment of this area, through 
including the potential formulation, the conservation of energy and momentum, the generalization 
for magnetic monopoles, as well as simplifying previously reported results through eliminating the 
need for the spacetime metric. 
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I. INTRODUCTION 



Maxwell's equations were first published in 186^ mathematically describing then recently 
discovered electromagnetic phenomena. His equations were written for three-space, requiring 
12 equations in 12 unknowns. These equations were later rewritten by Heaviside and Gibbs, 
in the newly developed formalism of dot and cross products, which reduced them to the four 
equations now seen in most modern textbooks^ and shown below in S.I. units 

V ■ E = ^, (Gauss' law); (1) 

V X B — \dtE = /ioJ, (Ampere's law); 
V X E + ^tB = 0, (Faraday's law); 

V ■ B = 0, (Gauss' law of magnetism). 



where E, B, J are conventional vector fields, with E the electric field strength V/m and B 
the magnetic field strength in units of Teslas and V is the three-vector gradient given by 

„ d d d 

V = ei— +62^ + 63 — , 2 

ox oy oz 



and Ci, 62, 63 are the three real Euclidean space orthonormal vectors, with Cj ■ Cj = 6ij, where 
Sij is the Kronecker-delta symbol. Maxwell's four equations given in Eq. ([T]) along with the 
Lorentz force law 

K = g(E + vxB), (3) 

completely summarize classical electrodynamics^. It can be seen that we have adopted the 
convention that a three-space vector v = f iCi + f 262 + "^363, is identified in bold, which is in 
agreement with ISO 80000-2 international standard for mathematical notation. 

Hestenes in 196fP, also refei^"^ in this journal, using the formalism of GA, further 
condense Maxwell's four coupled first order differential equations (DE's) into just a single 
first order DE 

^dt + v]F = cfioJ, (4) 



c 

where the field strength F = E + icB, with the sources held in J = cp — 3 . The symbols 
F, J, are called multivectors and as we shall later see can be visualized as containing a 



combination of scalar, vector, pseudovector and higher order terms. It might appear invahd, 
from the perspective of vector calculus, that in the source term J, we have subtracted the 
vector quantity J from the scalar quantity cp. However this is an unnecessary limitation, 
which severely limits the power of vectors. We can also see the benefit in relaxing this 
restriction, in that the field strength F can now hold both the electric and magnetic fields 
in a single variable. This is analogous to the formalism of complex numbers that contain 
both real and imaginary numbers as a single entity, which is a very successful formalism for 
two dimensions. See Appendix A to confirm the equivalence of Eq. Q with Maxwell's four 
equations. 

This comparison is investigated further in Table |T] where we compare GA, Gibbs vector 
calculus, and the tensor calculus formalisms. Tensor formalism requires two equations, four- 
dimensional spacetime and a metric in order to represent Maxwell's equations, compared 
with Eq. Q in three-space using GA, without the use of a metric. Tensor notation is briefly 
described in Appendix B in order to allow a direct comparison. Differential forms can also 
be used, but still require two separate equations to represent Maxwell's equations. 

Even though Gibbs was able to reduce Maxwell's twelve equations down to four, as men- 
tioned, his formalism for vectors had significant structural limitations. For example, the 
cross product only applies in three dimensional space, because in four dimensions there is 
an infinity of perpendicular vectors. However, probably most serious in terms of students 
learning physics, is that, conventional vectors do not integrate with established algebraic 
intuitions regarding basic operations. That is, there is no division operation, the cross prod- 
uct does not apply in two dimensions and one cannot freely add vectors to previously known 
algebraic elements (scalars), so that vector algebra becomes a monolithic structure unto 
itself. Hence the intuitive understanding of physics concepts, as well as general geometric 
understanding, which depends on the understanding of vectors, is significantly reduced. 

Historically, as vectors became more popular in physics and in various other fields, new 
scientific discoveries such as quantum mechanics and relativity meant that vector analysis 
needed to be supplemented by a basket of other mathematical techniques such as: tensors, 
spinors, matrix algebra, Hilbert spaces, differential forms etc. As noted irP, 'The result is a 
bewildering plethora of mathematical techniques which require much learning and teaching, 
which tend to fragment the subject and which embody wasteful overlaps and requirements 
of translation'. 
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EM theory 


Gibb's vector calculus 


Tensors 


OA 


Fields 


E,cB 




E + icB 


EM equations 


V • E = f, V X E + 9tB = 

V X B - ^dtE = /ioJ, V • B = 


Lip , ; 


UF = cfioJ 


Charge conservation 


V-J + lf = 


J% = Q 


□ • J = 


Energy in fields U 


ie(E2 + c2B2) 


e(F'>^F^ + iF^uFn 


Z 


Poynting vector S/c 


7r7(E X B) 


1 




Invariants 


c2B2 - E2 

^B E 

c 


iF'^^FaB 

2 ^h' 

UaB^sF'^^F^^ 




Action Ep — Ek 


pV-3-A 






Minkowski Force 


K = 7q(E + V X B) 


cKt" = qv^F'''' 


cK = —qvF 


Conservation energy 


^+V-S = -J-E 


cdpT'^^ = -g'^PFp^J^ 


cD • To = JF 


momentum S/c^ 


V- T =/9E + J X B 




ca-Ti = {JF)i 


Potential function A 


E = -Vy- ^, B = V X A 


Fap = Aa^p — Ap^a 


F = UA 


Lorenz gauge 


A = V-A+^f = 


A% = Q 


□ • A = 


Maxwell's Equations 






□^A = chqJ 


Stress energy tensor 


T (See below for definition) 


+ Ig'^^F^^F^" 


= \FnpF 



TABLE I: Comparison of mathematical formalisms, showing relative simplicity of mathematical 
expressions in GA. Note: The blanks in several rows in the column for GA is because the relevant 
equation is automatically part of the equation above. For example, there is a blank row below 
the F^ entry for GA, because the separate scalar and pseudoscalar components of this expression 
produce the two invariants in the vector column. The Maxwell stress energy tensor given by 
T-H- Tij = eo{EiEj — ^5ijE'^) + -^{BiBj — ^dijB"^). We have also the four-vector n defined as, 
no = Co and rii = Cjeo- 

Alternatively, the mathematical formalism of geometric algebra, developed by Clifford in 
1878, as convincingly argued by Hestene^, provides a single unified mathematical language 
for physics. Clifford's development of GA was based on work by Grassman on the exterior 
product and Hamilton's work on quaternions from 1843. With this formalism we can freely 
multiply and divide by vectors in a natural way as well as adding and subtracting diverse 
elements such as scalars and vectors. This is achieved primarily through combining the dot 
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product and cross product into a single new product called the geometric product, which 
for two vectors u and v is defined as 

uv = u- v + iuxv = u- v + uAv. (5) 

This simple device (which works in any number of dimensions) allows Maxwell's equations to 
be seamlessly melded together as shown in Eq. Q . Maxwell himself adopted quaternions in 
his follow up work in 1873, Treatise on Electricity and MagnetisrrF, writing his complete set 
of equations in this formalism. However Maxwell died in 1879 just after Clifford developed 
GA and so never had an opportunity to implement Clifford's approach, which subsumes 
quaternions into a complete algebraic structure. Maxwell's equations however, have stood 
the test of time and despite the fact that it is nearly 150 years since he first presented 
his equations, his basic premises and equations have remained valid. It has been pointed 
oulP^, for example, that Maxwell's equations assume a massless photon. If photons were 
not massless an extra current term would need to be added to the right hand side of his 
equations, but many precision experimental tests on the photon's mass have confirmed 
a mass approaching zero. Also Maxwell's third and fourth equations. Ampere's law and 
Gauss's law of magnetism, are a way of stating that magnetic monopoles do not exist and 
in fact despite diligent searches, their existence has not been verifiecP. Although this is 
consistent with special relativity, which asserts that B fields are created by relative motion, 
being demonstrated irP', by using electrostatics and relativity alone one can derive the 
Lorentz force law, Eq. ([s]). 

Ironically, also. Maxwell's equations developed forty years before Einstein's 1905 relativity 
paper, satisfy the relativity principle even though they were developed under the assumption 
of an aether. Hence, despite various attempts at modifications. Maxwell's original equations 
remain essentially intact. 

A. A description of the GA formalism 

The mathematical elements of geometric algebra (GA) are constructed from real scalars 
and unit vectors, with all other algebraic quantities being derived using the outer product. 
For example, for a three dimensional space we have the orthonormal basis vectors ei, 62 and 
63, satisfying ■ Cj = 6ij, with which we produce the bivectors Ci A 62, 62 A 63, Ci A 63 and the 
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trivector ei A 62 A 63. All algebraic elements can be freely added together to form a quantity 
called a multivector. 

Elements of the algebra form a linear space of multivectors over the real numbers. That 
is, if Ml and M2 are multivectors, then so is aiMi + 02^2 for any real numbers ai and 02- 

Multiplication between algebraic elements is defined to be the geometric product, which 
for two vectors u and v is shown in Eq. ([s]), where u ■ v vs the conventional symmetric dot 
product and m A f is the ant i- symmetric outer product related to the Gibb's cross product 
by M X = —iu A V OT u A V = iu X V, where i = 616263. 

The naturalness of this definition for the geometric product can be seen by using the 
distributive law to expand the product of two vectors 

(61M1 + 62M2 + e3U3){eiVi + e2V2 + e^v^) (6) 
= UiVi + U2V2 + U3VS + {U2V3 - V2U3)e2e3 + (uiv^ - U3Vi)eie3 + {uiV2 - viU2)eie2 
= u-v + uAv = u- v + mxv, 



using ef = 1 and 6j6j = —ejCi, which thus provides a definition of the dot and outer products 
in component form. This also shows how the basic distributive law for expanding brackets 
naturally splits into dot and outer products, thus integrating these with basic algebra. For 
a detailed calculation comparing the cross and the outer product using GA, please refer to 
Appendix C. 

For distinct basis vectors we find from Eq. ^ 

We can also see that i squares to minus one, that is = 616263616263 = 61626162 = — 1 
and commutes with all other elements and so has identical properties to the imaginary 
number j = The ubiquitous presence of a/— 1 in quantum mechanics, has been shown 

by Hestene^ to be due to the presence of spin. We thus give the unit imaginary 'i' its 
appropriate geometrical representation as the trivector i = 616263. 

The geometric product is in general not commutative but it is always associative, obeying 
a{bc) = {ab)c. We also find the very useful result that the geometric product of a vector v 
with itself, yields the scalar dot product, that is vv = v- v + zvxv = v- v = f^ + f| + f3. 



7 



Thus, from a pedagogical perspective, GA can more easily integrate with students' 
established intuition for elementary algebraic operations except that the concept of non- 
commutivity needs to be absorbed. However, it is actually a natural concept. For example 
the subtraction operation is non-commutative! That is5 — St^S — 5. A more physical 
example of non-commutivity is a series of rotations in three-dimensional space. 

In order to progress to four dimensional spacetime, we now add cq as the time axis, 
which anti-commutes with the space coordinates with Cg = 1. We then form a spacetime 
coordinate four-vector as 



where we denote a four-vector using the underline. This produces = (ct + x)eo(ct + x)eo = 
{ct + x)(ct — x) = c^t"^ — x^, with the correct invariant interval in spacetime. We could also 
define the product between four-vectors as xx^ in order to produce the same result. A four- 
vector is typically written v = [fo, "^^i, ^^2, ^^s] = [^^o^v], which is also the combination of a 
scalar with a vector, however in GA we provide the more general context for four-vectors as 
the two lowest orders of a multivector. 

We can see that the correct sign structure is produced through cq anti-commuting with 
each of the space axes 61,62,63. That is we do not need to add a metric to the basis, retaining 
a consistent Cq = = = e| = 1 with all vectors anti-commuting equally with each other. 
Thus we have achieved Minkowski's original vision of creating a four-dimensional spacetime 
with a complete symmetry between the space and time coordinated. 

So, summarizing our proposed notation for GA, following Hartld^^*^, we will use plain 
lower case letters for real scalars, (imaginary numbers are never used), plain upper case 
letters refer to general multivectors, unless indicated to be a vector by writing it in bold 
and with the unit three- vector written v, four- vectors represented with an underline and the 
unit trivector given by i. A general multivector can be written 



which shows in sequence, scalar, vector, bivector and trivector terms and can be used to 
represent many mathematical objects, such as scalars (a), complex numbers (a + ib), quater- 
nions [a + iw), vectors (polar) (v), four- vectors (a + v), pseudovectors (^w), anti-symmetric 
tensors (v + iw) with its dual (?v — w) and spinors (a + iw), with the four complex com- 
ponent Dirac spinor requiring the full multivector M. Refer to Appendix D for explicit 




(8) 



M = a + V + iw + ib 



(9) 



8 



mappings to the GA multivector representation, demonstrating how GA unifies a diverse 
range of mathematical formahsms. All products are the geometric product unless otherwise 
indicated. 

We now proceed to describe Maxwell's equations in GA, followed by the potential formal- 
ism that naturally follows and offers further interesting insights. We then describe the key 
equations used in solving problems in electromagnetism, such as the calculation of forces, 
and the conservation of energy and momentum. 

The use of the Lorentz transformations are then described, as well as the modifications 
required to Maxwell's equation for the existence of magnetic monopoles. 

II. MAXWELL'S EQUATIONS IN GEOMETRIC ALGEBRA (GA) 

We can re-arrange Eq. Q using the four-gradient □ notation as 

DF = cfXoJ, (10) 

where 

□ = i-^dt - V)eo, (11) 
F = E + icB = E^ei + E^es + E^e^ + ic{B^ei + ByCi + B^e^), 
J = (cp + J)eo, 



where i = 616263. The □ symbol represents the four-gradient, and we find = {^dt — 
V)6o(^9t - V)6o = i^dt - V){ldt + V) = jjdt - giving the d'Alembertian. It should 
be noted that there is a wide variation in notation within the literature when representing 
the d'Alembertian, such as A, □, or V; however following Feynman!^, we represent the 
d'Alembertian by and hence the four-gradient by □. 

For a moving charge distribution creating a current we have J = pv, where the four- 
velocity V = 7(c + v^ei + Vye2 + "^263)60 and we find = as expected. For free magnetic 
monopoles p™, we have the additional bivector and trivector current terms given by ip^v. 



The source current J, as defined in Eq. (11), is a four- vector, however we will not represent 
this with an underline in this case, because as was just demonstrated, with the inclusion 
of magnetic monopoles, J becomes a full multivector. With the addition of monopoles, 



the four-potential A also becomes a full multivector. From the four- velocity, we have the 
four- momentum defined as p = mv, giving = rn^c^. 

We also findD-J = — V)eo- (cp+ J)eo = |^ + V-J, so that setting the four-divergence 
to zero implies conservation of a quantity. 

Eq. (10) has been shown^ using Green's functions to have a solution 

where r' = |r — r'|, the distance from the field point r to the charge at r'. The prime denoting 
the value at the retarded time tr- Hence we can define a potential function 

A=^ [ -dT\ (13) 

giving F = UA. 



A. Potential formulation in GA 



When moving to the potential formulation, in order to simplify the equations, the Lorenz 
gauge is typically introduced (sometimes erroneously called the Lorentz gaugd^Sl). Heras et 
aP^ show that the Lorenz gauge is the only gauge that puts both the electromagnetic fields 
and the potentials on a causal basis with the sources. Also by multiplying both terms by gc, 
they take on the dimensions of energy, thus it can be seen that the Lorenz gauge expresses 
the conservation of energy. 

Looking at the Lorenz gauge in GA, we find it has the simple form of D-A = cV-A+^^ = 
0, as shown in Table [1} Hence the Lorenz gauge written in GA is equivalent to setting the 
scalar part of the geometric product to zero and to relax this restriction we simply use the 
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full geometric product, so that we then have simply 



F = UA (14) 

'1 



dt-Vj eo(V^ + cA)eo 
dt - V - cA) 



c 
1 

c 

c (^V ■ A + ^^t^^ + cV A A - - 9tA 
cA - - 9t A + icV X A 
cA + E + icB, 



where the Lorenz gauge A = V ■ A + Typically we set A = 0, however the presence of 



a non-zero Lorenz gauge could now be further investigated. Substituting into Eq. ( 10 ) we 
produce the second order D.E., 

W^A = c/ioJ, (15) 
where the four-vector A = [V^ + c^A^ei + Aye2 + ^263)160 = {V + cA)eo . This naturally 



splits into four copies of Poisson's equation 















A 





(16) 



which have known solutions, as shown in Eq. (13) 



The wave equation is generated from the case with no sources 

U^A = 0, (17) 

with each component satisfying the three dimensional wave equation {^d^ - V^)^ = 0. The 
general solution of which is 

(p{r,t) = ^{F{r-ct) + G{r + ct)), (18) 

where r = a/x^ + + z"^. The wave equation, including the polarization of light, can 
be further investigated using GA, as shown iii^. If we postulate that the photon is not 



completely massless, then we need to modify Eq. ( 15 ) as shown in Appendix E. 
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B. Forces in the fields 



The force on a charge moving in an electromagnetic field is given classically by the 
Lorentz force law, as shown in Eq. (|3]). We can write this in GA for a charge p moving with 
a four-velocity y_, in a field F as 

JF = pvF (19) 
= P7 (c + v) Co (E + iB) = p7 (c + v) (-E + iB) cq 
= 7 [-pv ■ E - cp (E - zv A B) + pv A E + ■ B] Co 
= -7 [pv • E + cp (E + V X B) + pv A E + ?v ■ B] cq, 



which we can see consists of scalar, vector, bivector and trivector terms inside the square 
bracket. The first term pv ■ E is representing the power density and the second term is the 
Lorentz force, hence we have produced the Minkowski four-force I^, so that we can write 

cK = -JF, (20) 



where we simply disregard higher order terms from Eq. (19) in order to produce the four- 
force i^. 

C. Conservation of energy - The work-energy theorem 

A key conceptual difference between Newtonian mechanics and electromagnetism is that 
instead of forces and energy conservation being connected to a localized object, these con- 
cepts now apply over a spread out field, and hence to apply the conservation of energy and 
momentum we now need to calculate a volume integral over the whole field, so that the total 
energy and momentum of both fields and matter is conserved. 

The work-energy theorem of electromagnetism, also called Poyntings' theorem, is nor- 
mally writterF' as 

-J-E = ^ + V-S, (21) 

where u the scalar potential energy and S = ^E x B = — ^E A B the Poynting vector, 
with each term having the units of power density. However, we can write this using the 
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four-gradient as □ ■ {cu + S)eo = eoi^dt + V) ■ {cu + S)eo = |^ + V ■ S, for the RHS, but also 
^FeoF = ^(E + zcB)eo(E + tcB) (22) 
= [-^(E2 + c^B^) + ^eo(EB-BE)]eo 

= [-^(^oE^ + -B')--ExB]eo 
= -{cu + S)eo 



and hence we can write Poynting's theorem as 

J-E = cn-To, (23) 

where 

To = ^eoFeoF (24) 

Also, we saw from the previous section that J ■ E is the scalar component of JF, therefore 
we can write Poynting's theorem in terms of F as 

cD-To = JF. (25) 

Therefore energy conservation for the free fields can be written 

□ ■ To = 0. (26) 

This is exactly parallel to charge p conservation with J = cp + 3 given by □ ■ J = 0. Hence 
with cTo = CM + S, we have energy u conservation given byn-To = ^ + V- S = 0. That 
is, if the energy u in the field falls, then this must be propagated away using the Poynting 
vector S. 

For example for a circular wire, shown in Fig. [l| of radius a carrying a current /, caused 
by a potential of V volts over its length L, we have E = and B = ^^e^. 
We find in cylindrical coordinates F = + i '^^°^2 ^4> and hence 

We notice that T has no time dependence hence we only require the divergence V ■ To and 
so the only term that will contribute the required scalar is 

Cr, (28) 



2ciia?L 
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FIG. 1: Variables defined for calculation of energy dissipated in a wire. 



which is the Poynting vector. If we calculate the divergence and sum over the volume 
we firstly find V ■ (re^) = 2 in cylindrical coordinates. Hence the power per unit volume 
is cD ■ To = and hence the total power for the volume na^L of the wire is VI, as 
expected. Instead of integrating the divergence over the volume we can use the divergence 
theorem, and sum over the bounding area. On the surface of the wire we have r = a and 



hence multiplying Eq. (28) by the the surface area A = 2TcaL we again find the power 



dissipated as VI. As a third alternative from Eq. (25), the power per unit volume is given 
by E ■ J = and hence the total power over the volume of the wire is also VI. 

The conservation of momentum using the electromagnetic stress-energy tensor is intro- 
duced in Appendix H. 



D. Conservation of energy using the four-potential 

It is stated irP^, that the essence of electromagnetism, can be expressed in the interaction 
energy J- A, where A is the four- vector potential. For example, for a loop of superconducting 
wire carrying a current, we have = and so ignoring the potential energy the total energy 
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is kinetic 



kin 



If A-Jdr (29) 

2 Jvol 



2 
1 

2 

i// A.dl 



JdaA • dl 

vol 



vol 



LP 



using ^ — LI, yielding the correct energy for an inductor. 
E. Lorentz transformations 

Simple experiments can confirm that moving a magnet past a wire will create an electric 
current or conversely a changing electric field will create a magnetic field. This can be viewed 
quite generally as a result of placing the observer into a frame of relative motion. 

When an observer moves past an electric field, it is found that the field parallel to the 
relative velocity is unaffected. Hence it is natural to split the field into two components 
parallel and perpendicular to the direction of motion, that is E = + Ej_. If the relative 
four-vector motion v between the frames of reference, is represented by a unit four-vector 
u = \veo = 7(1 + ^), then the transformed fields are found simply by appending this relative 
four-velocity onto the perpendicular components of the field, that is 

F' = E|| + E^ii. (30) 

For example, if we wish to calculate the electric and magnetic fields of a charge Q moving 
at a speed v/cin the ei direction relative to our instruments, then we can firstly look at the 
field from the simpler perspective of the frame at rest with respect to the charge. So, using 
Coulomb's law we have simply 

E = -^% (31) 
47reo 
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We then transform this into the moving observer frame, using Eq. (30), giving 



F' = E^ei + {Eye2 + ^.63)7(1 + -ei] 



Hence we have 



c 

V 

E^ei + j{Eye2 + E^e^) + -pi{Eye2 + -^^63)61 

V 

E^ei + j{Eye2 + E^e^) + -pi{-Eyeie2 + -^26361). 



(32) 



= E., = ^Ey, E'^ = ^E, 

B'^ = 0, cB'y = 7^E„ cS; = -i-^Ey, 



(33) 



which are the correct transformation laws for a pure electric fielcP, demonstrating how 
magnetic fields appear through relative motion. 

A similar approach works when translating coordinates of events between different ob- 
server frames, except that in this case we find that the perpendicular components are unaf- 
fected. Writing the four-coordinate a; = x_l + we then find the correct transformation is 
found by appending the relative velocity u onto the parallel components of the coordinates, 
that is 

X = + (34) 
For example for a boost in the ei direction of velocity v on coordinates a; = (ct + xei + ye2 + 



2:63)60, we find using Eq. (34) 



X 



{ct + X61) 6o7 (1 + ^ei j + (?/62 + 2:63) 60 



(ct + xci) 7(1 61 ) + 2/62 + 2:63 



eo 



[7 (ct - XV /c) + 7 (x - ft) 61 + j/62 + zes] 60, 



(35) 



which gives us the correct coordinate transformations 



, , XV \ 

ct' = 7(ct-— j. 



x' = 7 (x — vt) , y' = y, z' = z. 



(36) 
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More generally, if it is not convenient to extract parallel and perpendicular components 
from the field, then we can undertake the Lorentz transformations of the field or coordinates 
as a whole using a rotor. If we have a frame moving with velocity v — |v| in the v direction, 
then the field F transforms according to 

F' = RFR\ (37) 

where R — e~^"^/^ representing a boost in the v three-vector direction, where 9 is given by 
tanh 9 — v/c and hence 7 = cosh 9. We have also defined R^ — eoReo that corresponds to the 
Hermitian conjugate, where R reverses the order of products. The identical transformation 
law works to boost coordinates x = (ct+x)eo, being = RxR) . Thus the GA formalism has 

the advantage of using exactly the same transformation law for boosts, for both coordinates 
and fields. We can also do spatial rotations using the rotor S = e*"''^/^, when we are rotating 
about the axis of the unit vector u. 



III. SUMMARY OF KEY RELATIONS IN ELECTROMAGNETISM USING GA 

We can now summarize the relations required in GA to solve a typical set of electromag- 
netic problems. Given sources we have the field equations 

UF = C//0 J, (38) 

or in potential form 

W^A = cii^J (39) 

where F = E -|- icB = DA. With fields defined in hnear isotropic dielectrics we have 
F = c/xo(cD -|- iH), as seen in Appendix F. 
The conservation principles are 

□ •J = 0, Conservation of charge (40) 

□ • T = 0, Conservation of momentum and energy 

□ •A = 0, Conservation of energy (potential formulation) . 



We have the Minkowski four-force 

cK = -qvF. (41) 
17 



Through the use of the relativity principle we can work in a co-moving frame and convert 
four vectors simply by piggybacking the parallel component of the co-ordinate onto the 
relative velocity vector u 

X = Xy-u + (42) 

or alternatively, both fields and four- vectors can be transformed using the Lorentz transfor- 
mation 

a = RaR^ (43) 
where R = e~^"/^ and a can represent either a four-vector coordinate or the field F. 



IV. EXTENDED TOPICS 



We now present two specialized topics presented in the formalism of GA. An introduction 
to the Lagrangian formulation of electromagnetism using GA is also placed in Appendix G. 



A. Magnetic monopoles in GA 

One extension of Maxwell's equations that is commonly considered is the inclusion of 
magnetic monopoles. As mentioned, no free monopoles have yet been found, however a 
recent result in 2009 demonstrated an effective flow of magnetic monopoles in a spin ice 
systempSl^ though they cannot exist outside the material in a free form. There is an ongoing 
search for free monopoles at high energie^and so may yet be discovered. Hence the presence 
of monopoles should be allowed for by expanding Maxwell's equations. Fortunately these 



can be included very simply by expanding out the RHS of Eq. (10) to a full multivector 



that is, by including the bivector and trivector terms in the source J as shown, 

i^dt + V)(E + zcB) = c/io(cp - J - ^J'" + zp"^), (44) 

where p"* represents magnetic charge and J"^ is the current of magnetic charge. So we can 
see that we now have V ■ -B = /iop"^ and V x E + dtB = — /io J™' as required for the presence 
of monopoles. 



Thus the use of GA leads to obvious 'holes' in Eq. (10), which can very naturally be filled 



with the terms for magnetic monopoles in the multivector representing the source currents. 
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Also we can see that we can generalize no further because the multivector for J is now 
complete. 

Thus GA provides a natural framework for further investigation of the properties of 
monopoles. 

B. Lienard-Wiechert potentials 

For a single charge q in arbitrary motion with the four-vector distance between the charge 
and the field point r = (ct + r)eo, with four velocity u^. = (c + Vr)eo at the retarded time, 
we can write the four-vector potential as 

A(r,t) = -^^, (45) 
47reo r ■ 

where to place the charge on the past light cone, we specify = or c^t^ — = that 
implies |r| = ct. Therefore we find r ■ y_,,. = {ct + r)eo ■ (c + Vr)eo = c^t — r ■ = c|r | — r ■ v^. 



Thus 



1 ^lc + v^_ (46) 

47reo c\r\ — r ■ 



producing the well known Lienard-Wiechert potential^ 

Vir,t) = -r^^j^ , A(r,t) = ^y(r,t), (47) 

47reo c\r\ — r ■ 

where r is the vector from the observer to the moving charge at the retarded time tr given 
by |r I = |r — w(tr) \ = c{t — t^), where w is the position vector of the moving charge, moving 
at a velocity v^.. 

We could now calculate the E and B fields by differentiating these potentials, alternatively 
we can simply calculate the Coulomb field in a co-moving frame, and apply a Lorentz 
transformation to find the correct E and B fields, find the electric field 

E = ^7^4, (48) 
47reo 

where k is a scaling factor depending on the relative velocity and the angle of the velocity 
vector of the charge. 

So, for a charge in uniform motion, the electric field vector E as seen by the observer, there 
is an interesting quirk of nature, in that the electric field vector points towards where the 
charge is now instantaneously^^, indicated by f , not where the charge was when it generated 
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the field at the retarded time. This is a significant simplification for this special class of 
problems where charges and observers are in uniform relative motion, giving a convenient 
shortcut for students. 

V. DISCUSSION 

We have shown how the GA formalism provides a simple algebraic approach to solving 
several common situations encountered in the field of electromagnetism, while maintaining 
manifest Lorentz covariance in all our equations. The formalism can achieve this without 
the significant overhead of the identities for vector calculus or the index notation of tensors 
and its associated spacetime metric. 

We have also focused on a straightforward presentation for the calculation of forces and 
the equations for the conservation of energy and momentum, in order to allow an easy 
progression from experience gained in Newtonian mechanics. 

The potential formulation also neatly integrates into the presentation, which provides 
an alternate solution path for some problems. We adopt the Lorenz gauge that guarantees 
causality in all expressions. The Lorentz transformations take a particularly simple form in 
GA and allows easy use of the co-moving frame. 

The relative simplicity of GA is also illustrated visually by inspection of Table |T] and we 
would expect a much faster learning curve for students using the formalism of GA. It also 
should be noted that all results and equations are naturally in a full relativistic setting, and 
so no re-learning is required when relativistic effects are included. 

In higher dimensions than two, for example with spinors, it is typically difficult for stu- 
dents to visualize the complex numbers, but because in GA the imaginary has been 
replaced with the trivector i = 616263, results can immediately plotted on the three-space 
axes 61, 62, 63, allowing geometric intuition to be more easily developed. 

In order to introduce new students to the properties of electric and magnetic fields, it 
has been found that a series of basic experiments on the field properties are best undertaken 
and the relevant Maxwell equation introduced in sequence. Using GA, an overview can be 
given at each step, showing how each equation can be extracted from the master equation, 
Eq. ([4]), thus providing students with an overview, or a mind map within which new results 
can be placed. This would also be beneficial for past students to the field, who are familiar 
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with the physical experiments and can now be given a unified mathematical framework. 



VI. APPENDIX 



A. The equivalence of the GA form of Maxwell's equations 



From Eq. Q, Maxwell's equations in GA can be written out fully as 

+ + e2dy + e^d^ [E^ei + £^^62 + £"263 + c (53:6263 + ByC^iei + 5^6i62)](49) 



— - c/io ( J^ei + Jy62 + J263) 

Co 



This form has all dot products and cross products expanded in full, but formed into a single 
equation. This form reproduces Maxwell's equations without the use of a metric using just 
the three Euclidean space directions 61,62,63. 

Using the relations 6j6j = 5ij and the anti-commuting property CiCj = —ejCi, we can 
firstly equate scalar parts, to confirm that V ■ E = ^. Next, by collecting the 61,62 and 63 
terms we find 

^dtE^ei - eic{dyB^ - d^By) = -cfiojiJx (50) 
^dtEye2 - e2c{-dccBz + cd^B^) = -cfiol2Jy 
^dtE.e-s - e^c{d^By - cdyB^) = -cfio^fsJz, 

which when added give V x B — -^dtE = /ioJ as required by Ampere's law. We have the i 
trivector terms, giving 

ci{dA + dyBy + d,B,) = (51) 

giving V ■ B = 0. 

The 6162, 6163 and 6263 terms give 

e2e3dtB^, + e2e^{dyE,- d:,Ey) = (52) 
eseidtBy + eseii-d^E^ + d^E^) = 
eie2dtB^ + eie2{d:^Ey - dyE^) = 
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and adding these three equations we find V X E + ^tB = 0, thus Eq. Q reprod uces all four 
of Maxwell's equations. 

B. Introduction to tensor formalism 

Using tensors we can bring Maxwell's four equations down to two. We construct an 
anti-symmetric field tensor 





—Ex 



Ex Ey 

cB, 
-cB-, 



E, 

-CBy 

cB^ 



-E, cBy -cBx 



with the dual tensor G'^'^ = e'^^^^FaB and if we define the current density 



(53) 



Jti = Pov^ = 7Po(c, Vx, Vy, Vz) = {cp, Jx, Jy, Jz) = {cp, J), 



(54) 



then Maxwell's equations read 



pa0 



cpo 



(55) 



where the comma in front of a symbol implies differentiation with respect to this coordinate. 
That is we can define 

(56) 



and with the Lorenz gauge A^^ = we have 



(57) 



where the 4- vector potential A'^ = (V, cA^, cAy, cA^ 



If we have the proper velocity v,. 



dr 



7(c, Vx, Vy, Vz), where v is the velocity relative to 



an observer and 7 = 
force on a charge q 



then the 4-momentum is defined as j9 = mv. The Minkowski 



CLT CLT at 

Vectors transform as = A^a*^ and for the field tensor F^^'' = A'^A'^F^'^. 



(58) 
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We have the spacetime metric is given by 



1 














-1 














-1 














-1 



(59) 



C. Calculating the outer product 



The question often arises, of exactly how do I calculate the outer product? There are 
two ways to calculate the outer product. Firstly, if we are in three dimensions, then we can 
use the relation that u A v = m x v. Using the calculation method typically taught for the 
cross-product, using the determinant of the matrix 



ei 62 63 
Ml U2 Us 
Vi V2 V3 



we find 



U X V = {U2V3 - M3f 2)61 - {U1V3 - U3Vi)e2 + {uiV2 - ViU2)es. 

Hence, using i — 616263, we find 
u A V = iu X V 

= 6i6263[(li2t'3 " ^'2^i3)el - {U1V3 - U3Vi)e2 + {uiV2 - ViU2)e3] 

= {U2V3 - V2U3)e2e3 + {U1V3 - ^3^1)6163 + {uiV2 - ViU2)eie2. 



(60) 



(61) 



(62) 



More generally, we can find the outer product simply by expanding the components of the 
vectors as a geometric product and then deleting the scalar part, that is we find firstly for 
the geometric product 

uv = uiVi + U2V2 + U3V3 (63) 

+ {U2V3 - V2U3)e2e3 + {uiVs - 1^3^1)6163 + {uiV2 - ViU2)eie2, 
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using ef = 1 and CiCj = —ejCi, see Eq. (|6|. We know uv = u- v + uAv and that the dot 
product is a scalar, hence we simply remove the scalar part of the geometric product to give 
the outer product 

u A V = {U2V3 - V2Uz)e2e2, + {uiv^ - u^iVijCie^i + {uiV2 - viU2)eie2 (64) 

in agreement with result calculated using the cross product. However this approach, of 
expanding the brackets, works for vectors in any number of dimensions and so is more 
general as well as being more intuitive. 

Also, using the fact that the dot product for vectors is commuting and the outer product 
is anti-commuting we can write the specific relations 

a-b = -(ab + ba) (65) 
a A b = -(ab — ba). 



which enables us calculate these two products directly in terms of geometric products. 

Not just vectors, but general multivector products can be evaluated using the distributive 
law to expand the brackets and the anti-commutitive property of the orthonormal basis 
vectors CjCj = —ejCi and the orthonormal property = 1. As an example, the following 
product can be expanded as follows 

(63 + 26263 + 0(2 + 56162) (66) 
= 263 + 46263 + 2i + 5i - IO6163 - 563 
= -363 + 46263 - IO6163 + 7i. 



We can see how irreducible elements of the algebra (1, 6j, 6j6j, i), transform into each other 
and can be freely combined together when using GA. For example when solving Eq. (|4]), 
because the LHS, when expanded using the geometric product, can contain all eight possible 
algebraic elements, l,ei,eiej,i, then this single equation actually represents eight separate 
equations, which for a the field vector F = E + iB, will have six unknowns as the scalar 
and trivector terms are missing. This is similar to equating real and imaginary parts in a 
complex expression, except that in this case we have eight distinct elements to equate. 
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D. Mapping to the multivector 



The general three-space multivector in GA, shown in Eq. ([9]), can be used to represent 
a diverse number of mathematical objects. Firstly, we can obviously represent scalars (a), 
vectors (v) and complex numbers (a + ib), where i = 616263. However the quaternions i,j,k, 



defined by i^ = j" 



ijk 



-1, maps into the multivector iw, or component wise, 



i —J- iei,i — 7- — Z62,k —7- 263. Hence a general quaternion 

q = a + bi — ci + dk^a + ibei + ice2 + ide^ = a + iw. 
This is similar to the Pauli spinor mapping 



(67) 



-H- z/' = a + mi6i + ^262 + ^0363 = a + iw, 



(68) 



a + jas 
-a2 + jai 

also mapping to the even subalgebra of the multivector, which shows the close relationship 
between spinors and quaternions. A single qubit can also be represented as spinor in 

The electromagnetic field represented by the antisymmetric tensor F^'^, see Appendix B, 
maps as follows 

Ff"" ^F = E + icB, (69) 

with the dual tensor given in GA by G = iF. We also have other dual spaces in the 
multivector, such as for vectors \ -H- iw and for spinors a + iw ^ i{b — i\). 
The Dirac wave function maps to the full multivector as follows 

a + jas 
02 + jai 
-63 + jb 
-bi+ jb2_ 

Other entities can also be mapped to a multivector, such as pseudoscalars (ib) and pseu- 
dovectors (iw). Rotation matrices are also represented in a multivector as i? = a + iw, 
which rotates a vector about the w axis by an angle given by tan^^ = \w\/a according to 

y' = RyRI (71) 

Hence the great capability of the three-space multivector is demonstrated, being able to 
replace a large variety of mathematical structures and formalisms as shown and hence unify 
the mathematical language of physic^. 



^ ao + 6i6i + 62e2 + &363 + mi6i + m262+m363 + i6 = a + \ + iw + ib. (70) 
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E. Proca Equation in GA 



Maxwell's equations implicitly assume that the photon has zero mass and current pre- 
cision tests have concluded that its mass is indeed approaching zero, with an upper bound 
of < 6 X 10~^''eV^. However even a very small mass would be significant in that it 
would imply an additional longitudinal polarization direction for light as well as implying a 
variation from the inverse square law. 

Following ProcaF^, in order to correct Maxwell's equations, we add a mass term to Eq. ( 15 ) 
as follows 

□M = c/ioJ--^A (72) 

which in free space, J = 0, produces the Klein-Gordon equation. The photon's mass can be 
assumed to be a real mass, or an effective mass which can be imparted in plasmas. 



F. Quasi-static expansion of Maxwell's equations using GA 



The quasi-static expansion of Maxwell's equations was used by Chua in order justify 
the existence of a new circuit element he called the memristoi'^ZI. In the presence of linear 
isotropic dielectrics, Eq. Q can be written 

^dt + V]F' = J 



(73) 

where F' = cD + iH and J = cp — 3, where D = eE and B = fiH. We can define the vector 
fields, as a power series in a, for example using the electric field we have 



E = E 



Q=0 



+ a 



dE 

da 



+ 



a=0 



2 da^ 



a" d^E 



a=0 



Eq + aEi + a^Ea + ■ ■ ■ + a^'E^ + 



(74) 



where 



We can then write 



Ei 



1 d'^E 



k\ da'' 



0=0 



(^^dr + V) {Fq + aFi + a^F2 + ...) = (Jq + a Ji + + . . . ) 



(75) 



(76) 
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hence it is easy to see that the orders of the quasistatic expansion become 

VFo = Jo (77) 

drFi + VF2 = J2 



The process of solution is now very clear in GA. From the zeroeth order fields we calculate 
Fq, which is substituted into the first order fields to find Fi, and so on. 

G. Lagrangian formulation 

As an alternative to the Lorentz force law, we can derive the equations of motion from a 
Lagrangian, where L = T — U, with T representing the kinetic energy, and U the potential 
energy. The advantage of the Lagrangian approach is that just two scalar fields T and V 
are required, as opposed to the two vector fields, E and B as part of the Lorentz force law, 
shown in Eq. ([s]). The Lagrangian for electromagnetism can be written in GA 

L = Lfieid + L,^t = ^F' - J A, (78) 

with the interaction energy J ■ A = qV — q\ ■ A. Substituting F = OA we find 

L = ^{DAy - J A. (79) 

Lagrange's equations for a field in tensor notation are 

where refers to external forces. Thus in four-space we form four equations of motion, 
one for each coordinate. In GA, Lagrange's equations become 

For a particle free falling in a potential field we have Q = 0, thus substituting the Lagrangian 



Eq. (79) into Lagrange's equations Eq. (81), we find 



□ (— □a)-J = 0, (82) 
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and hence we recover 

= /ioJ, (83) 

very simply confirming the Lagrangian approach in GA. 

The Lagrangian of a relativistic charged particle of mass m, charge q and four velocity 
V, moving in an electromagnetic field with vector potential A we add a kinetic term to find 



L = —F^ -J A- -mc^. 
2/io 7 



(84) 



H. Stress-energy tensor 



Momentum conservation in the tensor formalism is normally written 



(85) 



where T^'' = eo^F^^'rjai^F^'^ + ^r]^^ F"^ Fap) is the electromagnetic stress-energy tensor. If 
written out in full we find 



rpfll/ 



U Sn,/C Sy/C Sz/C 

Sx/ C U Txx '-^xy Txz 

Syfc ^yx ^ ^yy ^y^ 

Sz/c Tzx ^ -^z; 



(86) 



where Tjj = e^EiEj + ^^BiBj. We note from Eq. (22), that the first row T^^ is given by 



^e^FcoF. In fact the next three rows can be shown to be 



Ti = FciCoF, 



where i G 1 ... 3. Hence the complete stress energy tensor is given bj^ 



T"" ^ FnF, 



(87) 



(88) 



where n= (1 + ei + 62 + 63)69 ^ or explicitly T^*^ = (Fn^F) ■ n^. Therefore we can now 



write Eq. (85) using GA as 



a-Ti = ki (89) 
for the three equations of momentum conservation, where the vector force components are 



given by Eq. (20) 
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